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ABSTRACT 

Coherent structures in the distribution of the Faraday rotation measure of extragalac- 
tic radio sources are isolated using wavelet transformation technique. A new algorithm 
of wavelet analysis for data points nonuniformly distributed on a sphere is developed 
and implemented. Signatures of the magnetic fields in the local (Orion) arm, the Sagit- 
tarius arm (and its extension, the Carina arm), the synchrotron Loop I and, possibly, 
the Perseus arm have been revealed using the RM catalogues of Simard-Normandin 
et al. (1981, 551 source) and Broten et al. (1988, 663 sources). Unlike earlier analyses 
of the RM sky, our approach has allowed us to assess the stability of the results with 
respect to modifications of the data sample. Only the aforementioned features remain 
stable under mild sample modifications. 

We consider separately low-latitude sources at \b\ < 10° and, using the model 
of electron density distribution of Cordes et al. (1991), we estimate magnetic field 
strength by comparing the model wavelet transform with that of the real data. Inde- 
pendent estimates of the mean magnetic field strength in the Orion arm using low- 
and high-latitude sources converge to 1.4 ± 0.3 fiG. Rotation measures of low-latitude 
sources provide a clear indication of a magnetic field reversal at a distance 0.6-1 kpc 
towards the Galactic centre. Our analysis has revealed for the first time the extension 
of the reversal in the Carina arm. Low-latitude sources from the catalogue of Broten 
et al. (1988) indicate a magneto-ionic structure in the direction of the Perseus arm 
with the magnetic field direction reversed with respect to that in the Orion arm. The 
extent of the region with reversed magnetic field near the Sun is 3 kpc or more in the 
azimuthal direction. The average pitch angle of magnetic field in the nearby spiral 
arms is 15°, and the mean field strength in the Sagittarius-Carina and Perseus arms 
is 1.7±0.3 /xG and 1.4± 1.2 /jG, respectively. The line-of-sight magnetic field in Loop I 
is estimated as 0.9 ± 0.3 /zG. 

We find firm evidence of a dominant even symmetry of the local mean magnetic 
field with respect to the galactic equator. Our results are compatible with a moderate 
large-scale north-south asymmetry, with magnetic field in the southern hemisphere 
being stronger in a region of at least 3 kpc in size. It cannot be excluded, however, 
that the asymmetry is local and results from vertical bending of magnetic lines in a 
region of about 400 pc in size, with the Sun being located close to the top of a magnetic 
loop whose magnetic field is 0.5 /xG stronger than the average field. 

Key words: magnetic fields - polarization - methods: data analysis - interstellar 
medium: magnetic fields - Galaxy: structure 



1 INTRODUCTION 

Faraday rotation measures of polarized extragalactic radio 
sources and pulsars are a direct and informative tracer of the 
large-scale magnetic field of the Milky Way (see the reviews 
of Spoelstra 1977, Zeldovich, Ruzmaikin & Sokoloff 1983, 
Vallee 1983a, Sofue, Fujimoto & Wielebinski 1986, Ruz- 
maikin, Shukurov & Sokoloff 1988, Wielebinski & Krause 



1993, Kronberg 1994, Beck et al. 1996, Zweibel & Heiles 
1997, Vallee 1997). Faraday rotation measure RM is a 
weighted integral of the longitudinal magnetic field along 
the line of sight, 



RM = K I n c H ■ ds , (1) 
Jo 
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where n c is the number density of thermal electrons, H is 
the magnetic field, L is the distance to the radio source and 
K = 0.81 rad m -2 cm 3 /iG -1 pc -1 (the positive direction of s 
is that from the source to the observer, so that field directed 
towards the observer produces positive RM). Since thermal 
electron density is relatively small in the intergalactic space, 
the main contribution to ([j]) for extragalactic sources comes 
from the Milky Way and from the radio source itself. As 
magnetic fields in different radio sources are uncorrelated, 
the distribution of RM in the sky for a representative sample 
of radio sources can be used to study the magnetic field in 
our Galaxy. 

For extragalactic radio sources, the effective value of L 
is determined by the angular coordinates of the source in 
the sky, whereas Faraday rotation measures of pulsars es- 
sentially depend on the pulsar position within the Galaxy. 
Therefore, Faraday rotation measures of pulsars and extra- 
galactic radio sources are most often considered separately, 
although one can envisage procedures of their joint analysis. 
In this paper we discuss the extragalactic RM sky and those 
features of the Galactic magnetic field that can be deduced 
from its analysis. We apply wavelet transform techniques, 
a method perfectly suited to isolate regular structures in a 
noisy signal. 



2 STUDIES OF THE RM SKY 

Early analyses of the configuration of the large-scale Galac- 
tic magnetic field using Faraday rotation measures of extra- 
galactic radio sources were performed by Morris & Berge 
(1964) (37 sources, not all of them extragalactic) and Gard- 
ner, Morris & Whiteoak (1969) (355 sources). These authors 
noted a strong asymmetry between the northern and south- 
ern Galactic hemispheres and concluded that the average 
magnetic field in the Solar neighbourhood is directed toward 
the Galactic longitude I = 70-80°. Morris & Berge (1964) 
claimed that the field has opposite directions above and be- 
low the Galactic midplane. Gardner et al. (1969) suggested 
that the apparent antisymmetry is due to a local anomaly 
at / < 60° associated with the radio Loop I. These au- 
thors also noted that the apparent direction of the field shifts 
to larger I with the increase of the latitude of the sou rces 
considered; we explain this behaviour below in Sect. |6.2[ 
Vallee & Kronberg (1973) measured rotation measures of 
252 sources (published by Vallee & Kronberg 1975) to esti- 
mate the direction of the field as I = 90° ± 10°, b = 5° ± 10° 
(with the region of Loop I excluded from the analysis) 
and concluded that the large-scale field is similarly directed 
above and below the Galactic plane. Ruzmaikin & Sokoloff 
(1977) considered 154 sources from the sample of Morris & 
Tabara (1973) and estimated the field direction as about 
I — 99° (with the Loop I region excluded) presuming that 
the field direction is the same above and below the mid- 
plane. Ruzmaikin, Sokoloff & Kovalenko (1978) considered 
the catalogues of Morris & Tabara (1973) and Vallee & Kro- 
nberg (1975) (with sources having |RM| > 100 radm -2 re- 
moved); these authors have determined the field direction 
in the southern Galactic hemisphere to be I — 91° ± 22° 
and stressed that it is not possible to determine reliably the 
global field direction in the northern hemisphere because of 
the local distortions due mainly to Loop I. 



Comprehensive catalogues of the Faraday rotation mea- 
sures of extragalactic radio sources were published by 
Tabara & Inoue (1980), Eichendorf & Reinhardt (1980), 
Simard-Normandin, Kronberg & Button (1981) and Broten, 
MacLeod & Vallee (1988), which contain 770, 346, 555 and 
674 sources, respectively. RM data were also obtained for 
several localized regions in the sky (see Vallee 1983a). The 
data used in this paper are discussed in detail in Sect. ^. 

Simard-Normandin & Kronberg (1979, 1980) used RMs 
for samples of 476 and 543 extragalactic radio sources, re- 
spectively. These authors averaged RMs of neighbouring 
sources within a 10° -15° radius of each source in order to 
reduce the noise. They also noted that the RM distribution 
in the southern Galactic hemisphere is significantly more 
ordered than that in the northern hemisphere. The largest 
localized feature is Loop I which produces a perturbation 
of about 30 radm -2 at b > +20° in the first and fourth 
Galactic quadrants. The pitch angle of the local large-scale 
magnetic field is about 15°, as estimated from the fact that 
the associated positive RMs of the sources at |6| > 5° have 
a peak at about I = 255° (instead of 270°) and RM changes 
sign at about I = 165° . These authors argued that a region of 
positive RM = 50-300 radm" 2 centred at I ~ 40°, b ~ +5° 
is due to a reversed magnetic field in the Sagittarius arm. 
These authors were unable to identify the reversal in the ex- 
tension of the Sagittarius arm to the fourth Galactic quad- 
rant (the Carina arm) using their methods. 

Andreassian (1980) considered 301 extragalactic sources 
with 251 of them taken from Vallee & Kronberg (1975) 
and the remaining 50 from Haves (1975) and excluded 
those having |RM| > 100 radm" 2 at 15° < \b\ < 75° and 
|RM| < 8 radm -2 , so that only 176 sources remained. He 
considered least-squares fits using a model for the local large- 
scale magnetic field to conclude that the field has opposite 
directions above and below the Galactic plane. Andreassian 
(1982) used 419 sources from the catalogue of Tabara & In- 
oue (1980) applying the same selection criteria (retaining 
the region of Loop I) to conclude that the fields revealed by 
the sources at b > 15° and b < 15° have opposite directions. 

Using the RM data published by Tabara & Inoue 
(1980), Inoue & Tabara (1981) estimated the field direc- 
tion as I = 100 ± 10° from 224 sources located at |6| > 30° 
and I ~ 80° from those at |6| < 30°. Sources with |RM| > 
100 radm -2 were excluded from their analysis. The original 
rotation measures were averaged over circles of a 15° radius 
centred on each source. These authors did not confirm the 
presence of a reversal of the regular magnetic field between 
the Orion and Sagittarius spiral armsQ and noted that the 
apparent reversal suggested by Simard-Normandin & Kron- 
berg (1980) may be due to Loop I. 

Sofue & Fujimoto (1983) also analyzed the sample of 
Tabara and Inoue (1980) consisting of 1510 extragalactic 
radio sources 770 of which have known RM, but applied less 
stringent selection criteria. They discarded an unspecified 
number of sources with |RM| > 300 radm -2 and with er- 
rors exceeding 10 radm -2 . The authors then used Gaussian- 
weighted mean values of RM for each source and studied 



* The feature called here the Orion arm is plausibly a relatively 
short spur between the more extended Sagittarius and Perseus 
arms, so that many authors call it the Orion spur. 
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features outside the local arm taking the averaging Gaus- 
sian with a full width at half-maximum of 20°. They found 
four alternating maxima and minima of RM along Galactic 
longitude at |6| < 20° and noted that they approximately co- 
incide with tangential directions to the inner Galactic spiral 
arms, and interpreted this as an evidence of the correspond- 
ing number of magnetic field reversals in the inner Galaxy. 

Attempts to determine the field direction in the next 
outer Galactic spiral arm, the Perseus arm, were made by 
Vallee (1983c) and Agafonov, Ruzmaikin & Sokoloff (1988) 
who used 45 sources located in this region of the sky. Al- 
though both papers used the same RM data, they came to 
opposite conclusions about the direction of the field in the 
Perseus arm: Vallee argues that it is parallel to that in the 
Orion arm, whereas Agafonov et al. conclude that the fields 
are almost antiparallel. 

Clegg et al. (1992) obtained and analyzed Faraday ro- 
tation measures of 56 low-latitude extragalactic sources at 
45° < I < 93°, \b\ < 5° and compared their longitudinal 
distribution with predictions of various magnetic field mod- 
els. They adjusted parameters of a circular-field model to 
argue that the field has a reversal between the Sagittarius 
and Orion arms producing a null in RM at / = 62° with RM 
mainly positive at I < 62° and negative at I > 62°. 

Han & Qiao (1994) used Faraday rotation measures of 
103 extragalactic sources at |fo| > 60° from the compilation 
of Broten et al. (1988) in an attempt to estimate the ver- 
tical magnetic field (sources with |RM| > 300radm~ 2 were 
excluded from their analysis). Their estimate is B z ~ 0.2- 
0.3 fiG with a direction from the south to the north Galac- 
tic pole. It is, however, unclear how one can reliably distin- 
guish the local north-south asymmetry associated, e.g., with 
Loop I from a global vertical magnetic field. 

A recent study of the regular magnetic field of the 
Galaxy using rotation measures of extragalactic radio 
sources was performed by Han et al. (1997). These authors 
used the catalogue of Broten et al. (1988), omitted sources 
with |RM| > 250radm -2 , averaged the rotation measures 
within circles of a 15° radius around each source and dis- 
carded those whose RM deviates strongly from the mean 
value in its circle. Altogether, 551 sources were considered. 
They noted that, in addition to the dominance of negative 
rotation measures at < / < 180° and positive RMs at 
180° < I < 360° resulting from the local large-scale mag- 
netic field, the distribution of RM in the first and fourth 
quadrants (-90° < / < and 270° < I < 360°) is apparently 
antisymmetric with respect to the Galactic equator (this is 
not confirmed by our results); they found a confirmation of 
the antisymmetry in the RM distribution for high-latitude 
pulsars (\b\ > 8°). These authors argued that this pattern 
(that includes the region of Loop I) results not from a local 
magnetic perturbation but rather from the odd symmetry 
of the global magnetic field in the thick Galactic disc within 
a few kiloparsecs from its centre. 

A general feature of the above analyses is that they 
provide rather divergent results concerning magnetic field 
directions in remote areas such as neighbouring spiral arms, 
although the estimates of the local regular magnetic field 
are more coherent. Furthermore, most authors who did not 
assume the even parity of the magnetic field with respect to 
the Galactic equator were unable to find convincing evidence 
of this symmetry in the observational data. These difficul- 



ties arise because it is unclear how sensitive the results are 
to the selection criteria, to the data sample chosen, and to 
the analysis technique applied, and because it is not easy 
to isolate the effects of the local perturbations associated 
mainly with Loop I. 

The above analyses were based either on fitting a cho- 
sen global model of the large-scale magnetic field to the ob- 
served distribution of Faraday rotation measures or merely 
on naked-eye adjustments of the model parameters. This ap- 
proach can be improved in the following respects. The data 
themselves are mostly sensitive to the magnetic field near 
the Sun because extragalactic sources with known RM are 
observed mainly at high galactic latitudes. Even nearby spi- 
ral arms (Sagittarius and Perseus) occupy in the sky regions 
as small as |&| < 5-10° in the directions tangent to the arms 
(assuming the scale height of 500 pc for the magneto-ionic 
layer). This difficulty cannot be completely removed by im- 
proving statistics and enriching the sample with low-latitude 
sources since it is difficult to distinguish remote large-scale 
features from nearby localized magnetic perturbations. (Pul- 
sar data are more favourable in this respect provided dis- 
tances to the pulsars are reliable.) Therefore, one can ex- 
pect significant difficulties in deciding between models of 
the magnetic field having similar properties near the Sun 
but differing in global features (e.g., field reversals along the 
galactic radius and axisymmetric versus bisymmetric geom- 
etry). Clearly, in this situation it is important to use well- 
substantiated, model-independent methods of quantitative 
statistical analysis rather than simple naked-eye fits. 

A natural way to avoid these difficulties is to study, as 
a first step, the magnetic field in a close vicinity of the Sun 
as this field is described by an irreducibly small number of 
independent parameters. Ruzmaikin et al. (1978) applied a 
regression method to obtain estimates of the strength and 
direction of the magnetic field near the Sun from Faraday 
rotation measures of extragalactic radio sources. However, 
they used first-generation RM catalogues which contain a 
rather small number of sources and many unreliable mea- 
surements. Ryadchenko & Sokoloff (1997, unpublished) ap- 
plied this method to the more recent catalogue of Simard- 
Normandin et al. (1981) to obtain basically the same results 
as Ruzmaikin et al. (1978). Appropriate statistical analyses 
were performed for pulsar data by Rand & Kulkarni (1989) 
and Rand & Lyne (1994) who fitted the data with models 
of local and global large-scale magnetic fields. 

To study the global structure of the Galactic magnetic 
field in an objective manner without prescribing any model 
of magnetic field beforehand, one may use methods similar 
to Fourier decomposition, that is expanding the observed 
spatial distribution of RM over a complete set of orthogo- 
nal basis functions. For spherical geometry, expansion over 
spherical harmonics seems to be natural. Two-dimensional 
spherical harmonic analysis using 65 extragalactic radio 
sources was performed by Seymour (1966, 1984) who con- 
cluded that the harmonics of the first four orders are statisti- 
cally significant, and that the dominant toroidal component 
of the Galactic magnetic field is symmetric with respect to 
the Galactic equator; the local direction of the magnetic field 
was determined to be from I = 112° to I — 292° (probably a 
misprint: the opposite direction was meant by the author). 
Seymour also concluded that there is a weaker antisymmet- 
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ric toroidal field increasing with distance from the Galactic 
equator, which the author identified with a halo field. 

Spherical harmonic analysis is not an optimal tool for 
studies of the global magnetic field of the Galaxy, again be- 
cause more or less distant parts of the Milky Way occupy 
small regions in the sky and would require the inclusion of 
high-order spherical harmonics into the fit. A generic prob- 
lem with Fourier expansion is that the basis functions (e.g., 
spherical harmonics) are spread over the whole sphere (i.e., 
they are not localized in space). Therefore, small-scale fea- 
tures are represented by a large number of phased Fourier 
modes. Statistical noise in the data and uneven distribution 
of the data points in the sky hamper reliable estimates of 
the required large number of the fit parameters. 

In this paper we apply wavelet analysis to the distri- 
bution of Faraday rotation measures of extragalactic radio 
sources in the sky. Our aim is to understand which partic- 
ular parameters of the Galactic magnetic field can be reli- 
ably extracted from the available data and then to estimate 
these parameters. Unlike previous analyses, we also verify 
the stability of our results with respect to modest modifica- 
tions of the sample in order to isolate those features in the 
RM sky which are less sensitive to the incomplete sampling. 
Our main goal here is not to present an interpretation of 
the observed RM distribution in terms of competing theo- 
ries of Galactic magnetic fields but rather to clarify which 
details in the RM sky are real and which can be artifacts 
of irregular data grids, incomplete sampling, diverse selec- 
tion criteria, etc. Results of the wavelet analysis of the RM 
sky and their interpretation in terms of interstellar magnetic 
fields are discussed in Sect. |^ and ^| 

Wavelet analysis (see, e.g., Farge 1992, Holschneider 
1995a) is a natural generalization of the Fourier analysis 
aimed at a consistent description of not only the scale but 
also of the spatial localization of structures. Unlike decom- 
positions in harmonic functions, wavelet analysis uses basis 
functions that are localized in both physical and wave num- 
ber spaces. This allows to represent localized structures in 
the data in terms of a tractably small number of parame- 
ters. Wavelet analysis is therefore a perfect tool for objec- 
tive studies of the RM sky. We discuss our implementation 
of wavelet analysis and in Sect. ^ and Appendix A. 

Wavelet analysis has been introduced in its present form 
by Grossmann & Morlet (1984) for analysis of seismic time 
series. Its physical applications are abundant, including the- 
ory of turbulence (Farge 1992), where this method was in- 
dependently introduced by Zimin (1981) (see also Frick & 
Zimin 1993), and now it is also used in astrophysics, e.g. 
for studies of large-scale structures in the Universe (Bijaoui, 
Slezak & Mars 1993; Malik & Subramanian 1997), cosmolog- 
ical velocity fields (Rauzy, Lachieze-Rey & Henriksen 1995) , 
solar (Nesme-Ribes et al. 1995; Lawrence, Cadavid & Ruz- 
maikin 1995; Halser et al. 1997; Frick et al. 1997b) and stellar 
activity (Soon, Frick & Baliunas 1999), and galactic mag- 
netic fields (Frick et al. 2000). 



3 DATA SELECTION 

The catalogue of Simard-Normandin et al. (1981) contains 
carefully selected values of RM for 555 unresolved extra- 
galactic radio sources. This catalogue is rather extensive and 



is notable for well-defined and consistent criteria of data se- 
lection. Most sources in the catalogue of Simard-Normandin 
et al. have |RM| < 200-300 radm -2 . It is generally believed 
that large values of |RM| of the order of, say, 500 radm -2 in 
such catalogues (based on observations with only moderate 
wavelength coverage) are hardly due to the contribution of 
the interstellar magnetic field but rather arise from magnetic 
fields in the sources, unreliable determinations of RM due to 
the mr ambiguity in polarization angles, nonlinear relation 
between the polarization angle and the wavelength squared, 
etc. (Ruzmaikin & Sokoloff 1979, Vallee 1980). 

The catalogue of Broten et al. (1988) (we used its up- 
date of 1991 kindly provided to us by Dr. J. P. Vallee) is 
based on similar methods and selection criteria; it contains 
data for 674 extragalactic radio sources. The catalogue con- 
tains a number of sources with very large |RM|, and many of 
them are at low Galactic latitudes. Together with the cat- 
alogue of Simard-Normandin et al. (1981), this catalogue, 
especially its low-latitude sources, provides the basis of our 
analysis. 

The catalogues of Tabara & Inoue (1980) and Eichen- 
dorf & Reinhardt (1980) contain a number of radio 
sources which are not present in the catalogues of Simard- 
Normandin et al. (1981) and Broten et al. (1988). Further 
results have been published for radio sources in selected re- 
gions in the sky: the Perseus arm (Vallee 1983b) , the Scutum 
arm (Vallee, Simard-Normandin & Bignell 1988) , the Mono- 
gem arc (Vallee, Broten & MacLeod 1984) and the Gum 
nebula (Vallee & Bignell 1983), and also for other sources 
at low Galactic latitudes (Clegg et al. 1992). Much of these 
data have been included into the catalogue of Broten et al. 
(1988). It is not a trivial matter to combine all the data into 
a larger joint catalogue because different authors use differ- 
ent selection criteria and procedures for the determination 
of RM from observations. In particular, the catalogues of 
Tabara & Inoue and Eichendorf & Reinhardt contain a sig- 
nificantly larger fraction of sources with very large values of 
RM than the catalogues of Simard-Normandin et al. and 
Broten et al. 

There are some contradictions between these cata- 
logues, so we have selected the data in the following man- 
ner. If the other catalogues disagree with that of Simard- 
Normandin et al. or Broten et al., we usually preferred the 
data from the latter two catalogues. Some exceptions were 
made when these two catalogues give a significantly larger 
value |RM| for a source than the other catalogues. We dis- 
carded such sources. 

The list of Clegg et al. contains many resolved sources 
with very large |RM|, and the authors could not find a 
unique value of RM for some of them. We neglected these 
data as well. After such a selection, only seven sources re- 
main in the list of Clegg et al. As it seems to be more im- 
portant to preserve the uniformity of the data set than to 
extend it insignificantly, we excluded also these seven esti- 
mates from our analysis. 

We also removed the sources which are not mentioned 
by Simard-Normandin et al. or Broten et al. and for which 
the other catalogues give different estimates of |RM| all ex- 
ceeding 100 radm -2 . However, if one of the catalogues gives 
an estimate lower than 100 radm -2 , we accepted it. 

The above procedure yields a joint catalogue of Faraday 
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rotation measures. The catalogue can be found at the URL 
http : / / plaxa. icmm .ru/~ lab4/RMjoint . dat . 



The total number of sources in the joint catalogue is 
841. Only 4 sources have been discarded from the catalogue 
of Simard-Normandin et al. and 41 source omitted from the 
catalogue of Broten et al. However, we appreciate that the 
joint catalogue still may contain inconsistent data. We also 
stress again the importance of the uniformity of the data set 
to be used. 

Our results are based on the catalogues of Simard- 
Normandin et al. (1981) (551 sources have been used) and 
Broten et al. (1988; 663 sources). Magnetic field parame- 
ters for the remote parts of the Galaxy have been derived 
from low-latitude sources represented mainly in the latter 
catalogue. We compared the results obtained from the cat- 
alogues of Simard-Normandin et al., Broten et al. and the 
joint catalogue to check the stability of our conclusions with 
respect to variations in the data set. We have also checked 
the stability of the results by removing sources with large 
|RM| from these catalogues. 

We consider only those magneto-ionic structures whose 
angular radius noticeably exceeds 20° (at least along Galac- 
tic longitude), the resolution limit of our technique (see Ap- 
pendix A). Analysis of structures at smaller angular scales 
using our approach would require a catalogue with a sig- 
nificantly more uniform and dense sky coverage. Therefore, 
we restrict our attention to the nearest spiral arms and the 
largest single radio structure, Loop I. These features are sta- 
ble with respect to modest modifications of the data. Struc- 
tures of a smaller angular size (more distant spiral arms and 
magnetic bubbles) can be revealed by specialized detailed 
analysis of selected zones (e.g., Vallee et al. 1988; Vallee 
1983a, 1993, 1997). 



4 WAVELETS: A SCALE FILTERING 
TECHNIQUE 

Wavelet analysis is based on a convolution of the data with 
a family of basis functions that depend on two parameters, 
the scale and location. This is a generalization of the Fourier 
transformation. In the latter, harmonic functions are used as 
a functional basis characterized by a single parameter, the 
frequency. In both wavelet and Fourier transforms, the basis 
functions have zero mean values. However, wavelet transfor- 
mations use oscillatory functions decaying at infinity, i.e. 
localized in space. This introduces an additional parameter, 
the position vector of the centre of a wavelet. In order to 
avoid excessive number of parameters, a family of wavelets 
is constructed from self-similar functions. The self-similarity 
can be implemented exactly in flat space (e.g., on a plane), 
but only approximately in curved space, e.g. on a sphere 
(Torresani 1995; Holschneider 1995b). 

The self-similarity distinguishes wavelet transformation 
from the windowed Fourier transformation, where the fre- 
quency, the width of the window and its position are in- 
dependent parameters. These basis functions are not self- 
similar and, in an isotropic case, contain three parameters 
(one of them, the window position, can be a vector), which 
complicates both interpretation of the results and the con- 
struction of the inverse transformation. Therefore, windowed 




Figure 1. The 'Mexican hat' wavelet in the (£, b)-plane (c = 2). 

Fourier transforms are usually applied only in the simplest 
one-dimensional case (e.g., to time series). 

We consider continuous wavelet transformation in two 
dimensions on a sphere of unit radius (the sky). We intro- 
duce angular coordinates / and b, the galactic longitude and 
latitude, respectively. The wavelet transform w(a, I, b) of the 
signal RM(/,fe) is defined as 



w(a, I, b) 



ip a , r (r')RM(r')(fr' 



(2) 



where ifj a ,r is the wavelet of a scale a centred at a position 
r = (I, b) on the sphere (r should not be confused with the 
radius vector) and the integration is over the unit sphere. 
The wavelet transform can be understood as the amplitude 
(positive or negative) of a structure of a scale a at a position 
(I, b). Wavelets must have zero mean value, 



ip a ,r(r') d 2 r' = 



(3) 



The choice of the wavelet depends on the goals of the 
analysis. For a local spectral analysis in time or space, 
wavelets with a good spectral resolution (i.e., well local- 
ized in Fourier space, or having many oscillations) should 
be used. To localize and analyze spatial structures repre- 
sented in the data, which is our goal here, a wavelet which is 
well localized in the physical space is preferable (i.e., rapidly 
decaying at infinity). 

We use a simple real isotropic wavelet known as the 
Mexican hat 

ip(x) = (c-x 2 )exp(-±x 2 ) , (4) 

where c is chosen to ensure that the mean value of ip(x) is 
zero; c = 1 in one dimension; the 'Mexican hat' in the plane, 
shown in Fig. [j], has c = 2. The wavelet family is generated 
from a function tp(x), called the analyzing wavelet, by its 
dilation and translation characterized by two parameters, 
the scale a and position r of a wavelet. A family of two- 
dimensional isotropic wavelets has the form 



i>a,r{r')=a 1 ip( y s{r,r')/a) 



(5) 



where r' = (l', b') is the current angular position in the sky 
and s is the angular distance between two points on unit 
sphere, r = (l,b) and r' — (l',b'), given by 



© 0000 RAS, MNRAS 000, 



6 P. Frick et al. 



2. O r 



1.0- 



30 60 

R, degrees 



90 



Figure 2. The integral wavelet spectra M(R) for the catalogues 
of Simard-Normandin et al. (1981) (solid), Broten et al. (1988) 
(dashed) and the joint one (dotted). The spectra are normalized 
by the maximum value of M(R) for the catalogue of Simard- 
Normandin et al. 



s(r, r ) = arccos[cos &cos b' cos (I — I 1 ) + sin b sin b] 



(6) 



We generalize the definition of a wavelet family by introduc- 
ing the dependence of c on a in Eq. (0) in order to implement 
an approximate self-similarity on a sphere (see Appendix A 
for details). Geometric factors in Eq. ^ make the wavelet 
anisotropic in terms of the angle variables (l,b). 

The traditional wavelet transformation should be fur- 
ther modified in our case because the data are given at ir- 
regularly distributed points in the sky, being neither a con- 
tinuous function nor given on a regular grid. Our implemen- 
tation, based on the gapped wavelet technique introduced by 
Frick et al. (1997a), is described in Appendix A. This mod- 
ification also affects c in Eq. (Q): now c becomes a function 
of the position of the wavelet centre, r. 

We present some of our results in the form of maps 
of the wavelet transform w(a, I, b) for a given scale a. The 
scale of a structure is characterized by the angular distance 
R from the centre of a wavelet to its first zero, i.e. ip\ x -R = 
and R(a) = c 1,/2 a - see Eq. (Q). R as a function of a is de- 
rived and discussed in Appendix A. The wavelet transforma- 
tion was performed for a wide range of R, but for practical 
reasons we varied R by steps of 3° (as determined by the 
typical separation of the data points) ; thus all the values of 
R quoted below have the accuracy of ±3° . We display only 
those wavelet transform maps which correspond to physi- 
cally distinguished values of R. The maps shown below are 
in fact scale-filtered RM distributions with a filtering win- 
dow centred at a certain scale R. 

Our aim is to isolate spatial structures in the RM dis- 
tribution in the sky and to estimate their parameters. For 
this purpose we consider the quantity 



M(R) = - 



w (R, r)d r . 



(7) 



In Fourier analysis, a similar quantity is known as the spec- 
tral energy density. The normalization factor in Eq. Q is 
chosen so that structures of the same amplitude but differ- 
ent scales would contribute equally to M (R) . As can be seen 
in Fig. ^, M(R) has a wide maximum at about R = 70° for 
all three catalogues used. The maximum in the spectrum is 
due to the largest structures in the RM distribution, that 
is the magnetic field in the Orion arm (see Sect. [5J). As in- 
dicated by a slower decrease of M with R, the catalogue of 




/, degrees 
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Figure 3. Contours of the wavelet transform ui(a, I, b) in rad m 2 
at the dominant scale R = 67° for the catalogue of Simard- 
Normandin et al. (1981). 



Broten et al. and the joint catalogue contain more structure 
at small and intermediate scales. 

We note that the wavelet used here is devised to obtain a 
good spatial resolution at the expense of a poorer resolution 
in the space of scales R (or the wave number space). There- 
fore, the imprints of underlying structures may be hardly 
visible in Fig. |^, the more so that M(R) is an integral pa- 
rameter which smears individual features even if they are 
strong. Such features are better visible in the wavelet trans- 
form maps in the (I, b) plane. 

The use of an isotropic wavelet does not preclude the 
detection of strongly elongated structures. In such cases the 
dominant scale R corresponds to the smaller dimension of 
the structure, and w at this scale has an elongated maxi- 
mum. The curvature of an elongated feature contributes to 
the wavelet transform at larger scales. 



5 STABLE FEATURES IN THE RM SKY 

The wavelet transforms presented in this section have been 
obtained from the sample of Simard-Normandin et al. (1981) 
unless stated otherwise. Figure ^ shows the wavelet trans- 
form at the scale R — 67° where M(R) is maximum in all 
three catalogues (within 3°). The main features of this map 
are two structures located at about / = 90° and 270°; they 
have opposite signs. These maxima are imprints of the large- 
scale magnetic field in the local (Orion) arm and Loop II 
(the Cetus Arc) (see also Sect. ^|). |RM| has slightly differ- 
ent peak values in the first and fourth quadrants and the 
zero-level contour is asymmetric in galactic longitude. The 
other two features visible in this map are a mild maximum 
of w at I — 30° and a minimum at I ~ 300° in the northern 
hemisphere. These are manifestations of Loop I (see below). 

It is notable that the dominant structures at I ~ 90° and 
270° are shifted from the Galactic equator to negative Galac- 
tic latitudes. The asymmetry between large-scale RM distri- 
butions in the northern and southern Galactic hemispheres 
is well known from earlier studies discussed in Sect. |^. 

In order to quantify the north-south asymmetry we con- 
sider for a moment the two hemispheres separately. In Fig. 4 
we show wavelet transforms and integral spectra for the two 
hemispheres. The spectra are maximum at very different 
scales, R = 35° in the north and R = 76° in the south. 
This indicates that the magnetic field in the northern hemi- 
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Figure 4. Wavelet transforms w(a, I, b) of the sample of Simard- 
Normandin et al. (1981) obtained separately for the northern 
and southern Galactic hemispheres at scales corresponding to the 
maxima in their respective energy spectra M(R). (a): contours of 
w (in radm -2 ) at R = 35° for the northern hemisphere (up- 
per panel) and at 76° for the southern hemisphere (lower panel), 
(b): M(R) for the southern (solid) and northern (dashed) hemi- 
spheres. Mo is the maximum of M(R) for the southern hemi- 
sphere. 



sphere is contaminated by distortions whose dominant an- 
gular radius is 35° . The large-scale features corresponding to 
the magnetic field in the Orion arm prevail in the southern 
hemisphere, but not in the northern one. 

There are two most conspicuous features in the north- 
ern hemisphere, which have no counterparts in the south 
and which are largely responsible for the asymmetry. One 
of them has the centre at (l,b) w (45°, 0), and the other is 
an elongated feature joining two local maxima at (l,b) ~ 
(335°, 30°) and (l,b) « (295°, 15°). These structures are in 
the region of Loop I (and Loop IV) and the Sagittarius arm. 
According to radio continuum observations at 820 MHz, 
Loop I occupies approximately the same region in the sky 
with its centre at (I, b) m (329°, 18°) and the angular di- 
ameter of about 116° (Berkhuijsen 1971). The remaining 
structures in the northern hemisphere visible in Fig. 4 have 
signs in agreement with those in the south, so that they can 
hardly cause a strong north-south asymmetry. 

In order to isolate the large-scale magnetic field in the 
Orion arm we first reduce the north-south asymmetry by 
omitting all the data in the region of Loop I. (Due to an 
only moderate scale resolution of our wavelet, the scale fil- 
tering, as in Fig. ^, is not sufficient to eliminate contamina- 
tion by smaller structures.) It is sufficient to omit points for 
/ from 280° to 50° at b > to match the dominant scales in 
the two hemispheres. The wavelet transform was then per- 
formed for the whole sphere assuming that the data in the 
above region are simply absent. Our gapped wavelet tech- 
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Figure 5. Filtered RM distributions: (a) at a scale R = 76° with 
the region of Loop I removed; (b), (c) and (d) at a scale R = 35° 
with the Loop I region retained but the large-scale contribution 
of panel (a) subtracted. The catalogue of Simard-Normandin et 
al. (1981) was used in (a) and (b), the joint one in (c) and that 
of Broten et al. (1988) in (d). Contours represent the value of 
w(R,l,b) in radm -2 . 



nique ensures a correct handling of such spatial distributions 
(see Appendix A) . 

A wavelet transform map for the dominant scale after 
removing the region of Loop I is shown in Fig. ^a; it is no- 
table that the dominant scale remains to be close to 76°, 
as for the southern hemisphere alone. The map in Fig. 
does not have a gap in the region where we have omitted the 
data because the wavelet scale is larger than the size of this 
region. It is instructive to compare Fig. |a with Fig. | which 
shows structures at about the same scale, but for the orig- 
inal data. The large-scale magnetic field in the Orion arm 
is responsible for the two extrema clearly visible in Fig. |^a. 
The spots are more symmetric with respect to the Galactic 
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equator than those in the original data. We use this map to 
study the large-scale magnetic field near the Sun in Sect. ^| 

In order to study structures at smaller scales, we sub- 
tracted from the original data the component of the scale 
R = 76° discussed above (i.e., obtained with the data in the 
region of Loop I omitted) . Then we wavelet-transformed the 
result on the whole sphere, including the region of Loop I. 
This procedure is meaningful provided the residual obtained 
after subtracting the dominant component can be proven to 
contain a significant signal above the noise level. We verify 
this by checking the stability of our results (including those 
at R ~ 30-40°) with respect to sample variations. Shown 
in Figs ^b,c,d are the wavelet transforms at R = 35° for 
the catalogue of Simard-Normandin et al. (1981), the joint 
one and that of Broten et al. (1988), respectively. The dis- 
tribution of the residual RM has a wide maximum of M(R) 
at R — 35-40° in all three catalogues. Apart from small- 
scale features, the filtered distributions of the residual RM 
at R = 35° shown in Figs|^b,c,d also reflect the global mag- 
netic field in the Orion arm; the corresponding structures 
are located slightly to the south of the Galactic equator. Es- 
pecially well pronounced and stable is that at I ~ 90°. Its 
counterpart in the fourth quadrant is weaker in Figs [|b,c 
(this is true also at R — 76°) and sensitive to the catalogue 
chosen: it can be more or less clearly seen at I » 250° in 
Figs ^3,d but is hardly visible in Fig. ^c. 

The structure of negative RM seen in Figs 5b,c,d at 
I w 90° , b < occurs in the Region A of Simard-Normandin 
& Kronberg (1979, 1980) identified with the Cetus arc (also 
known as the synchrotron Loop II, centred at I = 100°, b — 
—32.5° with a diameter of 91° - Berkhuisen et al. 1971, 
Berkhuisen 1971). The Cetus arc is known to produce a 
negative contribution to RM (Simard-Normandin & Kron- 
berg 1979, 1980). This makes it difficult to separate it from 
the contribution of the Orion arm that is also negative in 
that range of I. We have considered wavelet transforms at a 
range of scales 30° < R < 67° , but the difference in the scale 
between the two contributions (67° for the Orion arm and 
presumably about 45° for the Cetus arc) is not sufficient to 
separate them. Such a separation, based on an analysis of 
individual sources, is discussed by Vallee (1982). However, 
the shift of the wavelet extremum at I « 90° to negative lat- 
itudes can be, at least in part, due to the contribution of the 
Cetus arc. We note in this connection that, at scales about 
R = 35°, the negative extremum at I « 90° is significantly 
farther from the Galactic equator than the positive struc- 
ture at I ~ 250°. The RM structure due to the Cetus arc 
is plausibly better pronounced at smaller scales, but global 
methods similar to that employed here should await new RM 
catalogues to isolate the structure reliably and objectively. 

The positive extremum visible in Fig. gp at I fn 250° 
arises in the region which is close to the Gum nebula centred 
at I = 260°, 6 = (Vallee & Bignell 1983, Vallee 1993). The 
data used here have not allowed us to reliably identify this 
feature as its radius, ~ 20° is close to the resolution of our 
analysis. This feature is largely due to a few strong sources in 
the catalogue of Broten et al. located in this area. A possible 
signature of the Gum nebula can be seen as an enhancement 
in the magnitude of the positive extrememum of the wavelet 
transform to 35radm -2 at scales R < 40° (as compared to 
20 rad m -2 at larger scales) , indicating a narrow intense peak 



on top of a wider maximum. The position of this extremum 
at I ~ 250°, b ~ —10° changes only slightly as R varies. 

There are more details of a scale R — 35° located in 
the northern hemisphere. Clearly seen in Fig. ^b are a pos- 
itive extremum centred at (I, b) = (50°, 21°) and a nega- 
tive one at (l,b) = (289°, 23°). The position and intensity 
of the former structure are different for the catalogue of 
Simard-Normandin et al. and the joint one where it shifts to 
(l,b) = (110°, 41°) and its intensity changes by 25 per cent. 
The position of the structure centred at (l,b) = (289°, 23°) 
remains relatively stable in all three catalogues. This feature 
is apparently connected with Loop I and, possibly, Loop IV 
which can be seen at about the same position in radio contin- 
uum surveys (Berkhuijsen, Haslam & Salter 1971; Berklmi- 
jsen 1971; Salter 1983). 

The pattern of the signs of RM seen in Figs [|b,c in 
the general direction to the Galactic centre (with RM sign 
changing as — | h when moving clockwise from the quad- 
rant 360° > I > 270°, b > 0), was interpreted by Han et al. 
(1997) as indicating an odd (dipolar), axisymmetric mag- 
netic field in the central region of the Milky Way. However, 
the wavelet transform maps presented here indicate that this 
pattern is just a part of a more general RM distribution that 
extends wider along Galactic longitude and is more plausi- 
bly due to local magnetic field distortions. The significant 
north-south asymmetry is also better compatible with the 
local nature of the RM features at these scales. 

The available data on Faraday rotation measures of ex- 
tragalactic radio sources cannot provide stable and selection- 
independent results for magnetic structures at angular radii 
below 35° when a global analysis of the kind performed here 
is applied. This scale is equal to about seven times the typ- 
ical separation of data points in the catalogue of Simard- 
Normandin et al. (1981). As we show in Appendix A, 20° 
is the resolution limit for this catalogue (and for that of 
Broten et al. 1988) associated with the uneven distribution 
of the sources in the sky and their limited number. To check 
the robustness of our results, we also considered the result 
of removing sources with the largest |RM| in the catalogue 
of Simard-Normandin et al.; the structures at R ~ 30-40° 
remained stable (so that a structure does not shift by more 
than its radius) even when 19 sources with |RM| > 200 have 
been removed (of the total of 551 sources). Sources with 
large positive RM occur mostly near I ~ 220° and those with 
large negative RM are concentrated near I ~ 90°, therefore 
they are likely to trace real structures of a small angular size. 
The signs of RM in these two groups of sources are consis- 
tent with the direction of the local large-scale magnetic field, 
although widely distributed sources with moderate RM also 
provide a reliable, clear evidence of this field. The results 
obtained from the joint catalogue are less stable even at 
large scales, mainly because its data points are distributed 
less evenly in the sky. (We emphasize that the omission of 
sources with large |RM| was only used to assess the relia- 
bility of our results; our conclusions reported below are all 
based on samples that do include sources with large |RM|.) 

The resolution limit of the sample somewhat hampers 
analysis of magnetic field features located far from the Sun 
because they span a narrow range in b. Still, one can study 
some such details, provided they are extended along lon- 
gitude, by using only low-latitude sources and performing 
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a one-dimensional wavelet t rans formation in longitude. We 
discuss such results in Sect. 6.2. 



Since the scale resolution of the wavelet family used 
here is only moderate, we have restrained ourselves from 
discussing structures at intermediate scales, but instead con- 
sider only two widely separated scales, 76° and 35° both of 
which are distinguished by maxima in the wavelet spectral 
energy density. 

To conclude this section, we present in Table 1 pa- 
rameters of the stable structures discussed in Sect. ^| where 
magnetic field strengths also are estimated. The mean val- 
ues of the parameters for the Orion arm and Loop I have 
been obtained using the catalogue of Simard-Normandin et 
al. (1981). Their errors (except for magnetic field strength) 
were obtained as a maximum deviation, from the above 
mean value, of the values obtained from the joint cata- 
logue and pairwise combinations of the catalogue of Simard- 
Normandin et al. (1981) with those of Tabara & Inoue (1980) 
and Eichendorf & Reinhardt (1980). When estimating the 
errors, we omitted sources with |6| < 10° from the latter two 
catalogues and the joint catalogue. Furthermore, we added 
squares of the above maximum deviations and of the average 
angular separation of the measured points in the catalogue 
of Simard-Normandin et al., equal to 4°, to obtain the fi- 
nal errors of lo and bo. In a similar manner, the error of 
W given includes the uncertainty of the normalization con- 
stant c = h — K, (see Eqs (AS) and (A3)), which is equal to 
2 rad m -2 for the catalogue of Simard-Normandin et al. Indi- 
vidual features are further discussed in Sect. ^. The feature 
labelled Loop I in Table 1 occupies only a part of Loop I as 
described by Berkhuijsen (1971); in particular, the brightest 
part of the radio Loop I, the North Polar Spur, is slightly 
shifted from the above feature of the RM distribution. 

The results obtained from the catalogue of Broten et al. 
agree with those from the data of Simard-Normandin et al. 
to within 1/3 of the errors quoted in Table 1. 



6 THE LARGE-SCALE MAGNETIC FIELD 

In this section we discuss the magnetic field in the Galactic 
spiral arms and in Loop I considering separately high- and 
low-latitude sources. We estimate ma gnetic field strength in 
the Orion arm and Loop I in Sect. S.l and then consider low- 
latitude sources in Sect. 



6.2 to study magneto-ionic struc- 



tures in the Sagittarius, Carina and Perseus arms. 



6.1 The Orion arm and Loop I 

In order to estimate the strength of the large-scale mag- 
netic field in the Orion arm we compared the wavelet trans- 
form for RM in the whole range of b from the catalogue 
of Simard-Normandin et al. (1981) with those for model 
magnetic field distributions and adjusted B to obtain, at 
R = 76° , structures of the same intensity w and localization 
(lo, bo) as in the observed RM map. It is important to stress 
that the model wavelet transforms were always calculated 
on the real data grid (551 points positioned as in the cat- 
alogue of Simard-Normandin et al.). We have verified that 
the results change insignificantly when sources within the 
belts \b\ < 5°, 10° and 15° are omitted. 

The model magnetic field was purely azimuthal with 



strength independent of both galactocentric radius and the 
vertical coordinate. The half-thickness of the magnetic field 
distribution was adopted as h = 1000 pc (twice the scale 
height of the warm interstellar gas) and its radius as R s = 
15 kpc. We considered models with uniform thermal electron 
density and also with the electron density distribution of 
Cordes et al. (1991, see also Taylor & Cordes 1993), 



1 cm -3 



"2 



1 cm -3 



0.025 exp 



= 0.2 exp 



1 kpc V 20 kpc 



0.15 kpc 



4 kpc 



2 kpc 



(8) 
(9) 

(10) 



The radius of the solar orbit was adopted as Rq = 8.5 kpc. 
The term n\ describes a smooth disc, and ri2 represents a 
dense annulus at r = 4 kpc. Note that Heiles, Reach & Koo 
(1996) argue that this annulus does not exist. We also tried 
h = 2000 pc and i? g = 20 kpc for the size of the magnetic 
disc to check that our results are affected only insignifi- 
cantly in models with n c given by Eqs ([jj|)-(flb]) . We con- 
sider a purely azimuthal magnetic field in these fits because 
the quality of the data does not warrant more complicated 
models with non-zero radi al fiel d. The data for low-latitude 
sources discussed in Sect. |6.2.2 do not admit such a simple 
model, and so allow us to estimate reliably the magnetic 
pitch angle. 

We calculated model RM for all 551 points and then 
performed their wavelet transform. For a uniform disc model 
the results are as follows. The maxima of |w| occur at the 
same scale R = 76° as for the real data. The value of mag- 
netic field was found by adjusting the value of the model 
wavelet transform w at the extrema to those for the cat- 
alogue of Simard-Normandin et al.; B m 1.2 ^tG fits both 
the extrema. The main deficiency of this model is that the 
extrema of the wavelet transform are displaced from the ob- 
served structures towards I — by about 20° . 

Therefore, we further tried a nonuniform electron den- 
sity distribution retaining only m, Eq. (^|). Then the po- 
sitions of the extrema agree much better and, estimating 
B as above, we obtain 1.4 and 1.5 /iG for the positive and 
negative extremum, respectively; these values are cited in 
Table 1. The difference of the values of B obtained in vari- 
ous models gives the uncertainty of B as shown in the table. 

The annulus described by ri2 in Eq. ([[(^) is projected 
onto a small region at \b\ < 10° and so it does not affect 
rotation measures at large scales. 

The average line-of-sight magnetic field in Loop I can 
be estimated as 0.9 ± 0.3 jj,G from the contribution of Loop I 
to the Faraday rotation measure given in Table 1, \W\ = 
28 ± 8radm -2 , using the ionized gas density in the shell of 
0.4 cm" 3 and pathlength of 100 pc (Heiles et al. 1980, Salter 
1983). The amplitude of the negative extremum identified 
with Loop I, W = —28 ± 8radm -2 is close to the esti- 
mate of Faraday rotation measure in the North Polar Spur, 
~ +26 rad m -2 (Heiles et al. 1980), but the signs are dif- 
ferent. The latter estimate refers to a dense shell of Loop I 
at (l,b) — (32°, 42°), whereas the former to its part on the 
other side of the shell centre. The ensuing geometry of mag- 
netic field in the shell of Loop I is compatible with interstel- 
lar horizontal magnetic field bent into an arc and compressed 
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Table 1. The dominant structures in the RM distribution 



Name 


I,. 


00 


J? 


viz 


L> 


7*rcv 




(°) 


(°) 


(°) 


(rad m — 2 ) 


(AtG) 


(kpc) 


Sagittarius arm 


33 ±5 




35 


30 ±3 


1.7 ±0.3 




Carina arm 


330 ±5 




35 


-30 ±8 


1.6 ± 0.3 


7.9 














Orion arm 


98 ± 18 


-13 ± 17 


76 


-27 ±4 


1.5 ± 0.3 




Orion arm 


276 ± 17 


-17±9 


76 


24 ±8 


1.4 ±0.3 


10.5 














Perseus arm 


130 ± 30 




< 30 


20 ± 10 


1.5 ± 1.2 




Perseus arm 


190 ±5 




< 30 


-30 ± 10 


1.2 ± 0.7 




Loop I 


289 ± 16 


23 ± 10 


41 


-28 ±8 


0.9 ± 0.3 





Notes: lo and bo are the Galactic coordinates of the centre of a structure, R is 
its angular radius (scale), and W = w(R,lo,bo) is its amplitude; B is the total 
strength of the mean magnetic field, except for Loop I where the line-of-sight 
component is given; r rov is the galactocentric radius of a magnetic field reversal 
between the corresponding arms at I = assuming the radius of the Solar orbit of 
8.5 kpc. Negative (positive) values of W indicate that the field is directed toward 
(away from) the observer. Results for the Sagittarius arm (and its extension to 
360° > / > 270°, the Carina arm) and the Perseus arm have been obtained from 
low-latitude sources, |6| < 10°. 



by expanding shell, so that the line-of-sight magnetic field 
has opposite directions on the opposite sides of the shell 
centre. The North Polar Spur itself has a scale smaller than 
those discussed in this paper, but the structure of positive 
RM at (I, b) = (50°, 21°) in Fig. 5b can be related to it. Our 
estimate of magnetic field in Loop I agrees with the results 
of Troland & Heiles (1982) who argue in favour of a weak 
magnetic field enhancement in the shell of Loop I. 

The RM features attributed to the reversed magnetic 
field of the Sagittarius arm occur at about the same Galac- 
tic longitudes (and have the same signs) as Loop I and its 
most prominent part, the North Polar Spur. The connection 
of these features with a remote region is confirmed by the 
studies of pulsar Faraday rotation measures mentioned in 
Sect. ^| It is also noteworthy that the North Polar Spur is 
not well pronounced at low latitudes b < 8° (Berkhuijsen et 
al. 1971) where the structure attributed to the Sagittarius 
arm is well visible. 

To conclude this section, we give a crude estimate of the 
size of the region sampled by the extragalactic sources. Even 
though all sources have been taken into account to estimate 
B, the main contribution to RM (and the wavelet transform) 
is due to a region within a distance Lb = W/(Kn a B cos a), 
where n c — 0.03 cm" 3 is the typical value of electron den- 
sity and a is the average angle between the magnetic field 
direction (assumed to be from I — 270° to I = 90°) and 
the direction to a source. This yields 2Lb — 3 kpc for the 
size (along the field direction) of the region sampled. The 
low-latitude sources sample a more extended region. 



6.2 Low-latitude sources and magnetic field 
reversals 

As discussed in Sect. ^ the resolution limit of the cata- 
logues used here is about 20°. Therefore, features at |&| 
10° cannot be studied using an isotropic two-dimensional 
wavelet. Hence, we considered low-latitude sources sepa- 
rately to study remote structures located outside the Orion 
arm whose radius is larger than 20° in I, but smaller than 
10° in b. 



The large-scale magnetic field in the Orion arm pro- 
duces rather different structures in the RM distributions at 
low and high Galactic latitudes. Sources at latitudes \b\ > b, 
sample magnetic field in a region around the Sun whose ra- 
dius is L < h/ sin 6*. Under a crude assumption of a purely 
azimuthal field, the curvature of the field lines of the large- 
scale magnetic field can be neglected provided L < Rq, or 
6| > 7° for h = 1 kpc and Rq = 8.5 kpc. For sources at lower 
Galactic latitudes, \b\ < the curvature of the magnetic 
lines cannot be neglected. We adopt the borderline latitude 
as 6. = 10°. 

Hence, maximum values of |RM| produced in the Orion 
arm should be observed at I ~ 90° and 270° for sources at 
\b\ > 10° even for curved magnetic lines, but for low- latitude 
sources the largest |RM| are shifted towards I — 0. There- 
fore, high- and low-latitude sources should be considered 
separately. 



6.2.1 The catalogue of Simard-Normandin et al. 

The catalogue of Simard-Normandin et al. (1981) contains 
48 sources with \b\ < 10° which are shown in Fig. The 
scatter of the points is explained, in part, by scatter in the 
latitudes of the sources. We show in Figs qb-d the values 
of RM for the same 48 positions, but calculated from three 
model distribution of magnetic field B and electron density 
n e described in what follows. In order to provide a measure 
of the variations of RM with b, we present in these figures 
the model values of RM at both 6 = (solid) and |6| — 
10° (dashed). Note that the maximum values of |RMj may 
occur away from b = if the field has a reversal because the 
maximum value can be attained on the line of sight with 
\b\ 7^ that passes through the Orion arm but not through 
a more distant region with reversed magnetic field. 

The random Galactic magnetic field leads to a random 
walk of RM, and the associated random scatter of RM is 
especially strong for low-latitude sources where the number 
of turbulent cells along the line of sight is large. The data 
are further contaminated by the intrinsic rotation measures 
of the radio sources. To filter out the random fluctuations 
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Figure 6. The longitudinal distributions of RM for low-latitude 
sources, |6| < 10°, (a) in the catalogue of Simard-Normandin 
et al. (1981) and for the following four models: (b) a uniform 
azimuthal magnetic field and n c = ni, see Eqs (g)— (E(j), (c) a 
uniform azimuthal field and n c = n\ + n2, and (d) an azimuthal 
field with a reversal at r = 7.9 kpc in a disc with n e = m. The 
model RM distribution at b = is shown with solid line and 
that at b = 10° is shown dashed; dots show the values of RM, 
as observed (panel a) and calculated for the true positions of the 
sources (panels b— d) using the model magnetic field. 



that dominate at small scales, we compare with each other 
the wavelet transforms of the observational data and the 
model distributions, shown in Fig. ^, rather than the original 
distributions of RM along I shown in Fig. ^. 

We performed one-dimensional wavelet transformation 
(along I) for the data in Fig. ^| and show in Fig. the wavelet 
transforms in the (I, R) plane. The one-dimensional 'Mexi- 
can hat' was used as the analyzing wavelet. Unlike maps in 
the (I, b) plane discussed above, Fig. Inconveniently displays 
all the scales simultaneously. Figure[7k shows the wavelet 




180 



/, degrees 



Figure 7. One-dimensional wavelet transforms of the distribu- 
tions shown in the respective panels of Fig. ^ shown in the (I, R) 
plane with R the scale of the structure. Contour labels are given 
in rad m — 2 . 



transform for the data of Simard-Normandin et al., Fig. ga. 
The following three model distributions are presented in 
Figs ^|b-d and their wavelet transforms are shown in the cor- 
responding panels of Fig. [jj. The first model (Figs [jb and ^b) 
is a disc of a radius 15 kpc with a uniform azimuthal mag- 
netic field of a strength B (an adjustable parameter) and 
the electron density given by Eq. . The positions of max- 
imum |RM| in this model are strongly shifted towards I — 
and occur at I « 40° and 320° . The asymmetry between the 
left- and right-hand sides of these and other similar plots 
is entirely due to asymmetry in the data point distribution 
where more sources at < I < 180° happen to probe regions 
with large |RM|. We evaluated B by minimizing the r.m.s. 
deviation of w from the observed values of Fig. Ma to obtain 
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B — 0.8 fiG. These results change insignificantly if we take 
n e = 0.03 cm -3 independently of radius. 

The second model (Figs ^3 and ^jp) has a uniform az- 
imuthal magnetic field, but the electron density contains 
both terms of Eq. (^). This model results in very strong 
extrema of RM located at I w 25° and 335°. Only a few 
sources in the sample are projected on the annulus ni be- 
cause it occupies a narrow range \b\ Ss 1° at I m 25° and 
335°. 

The above models produce only two clearly visible max- 
ima of RM, whereas the observed picture, Fig. Wa, has a more 
complicated structure with four local extrema. Therefore, we 
had to consider a more complicated model that uses again 
n e = m, but contains a reversal of azimuthal magnetic field 
at a certain galactocentric radius r = r lev (the results shown 
in the figures correspond to r rcv = 7.9 kpc). The wavelet 
transform in Fig. 0d contains an additional pair of struc- 
tures produced by the reversed magnetic field. The model 
of Fig. 0d clearly provides a fair fit to the data of Fig. 0a. 
We adjusted both B (assumed to be uniform apart from the 
reversal) and r rov to reach the best agreement of the model 
and the 'observed' values of w by minimizing the r.m.s. de- 
viation. A minimum deviation is achieved for r rcv = 7.5 kpc 
and B — 1.4 fiG. However, the best overall agreement in 
terms of the shapes and positions of all the structures oc- 
curs for r rcv = 7.9 kpc and B = 1.3 fiG. Both estimates of B 
in the Orion arm agree within errors with that from RMs of 
high-latitude sources. Magnetic fi eld in the Sagittarius and 
Carina arms is estimated in Sect. |6.2.2 . 

A notable feature of Fig. 0d is that the field reversal in 
the inner Galaxy can be clearly seen in the wavelet trans- 
forms of both the original and model data in both first and 
fourth Galactic quadrants. There are only a few sources with 
negative RM in the Carina arm region (the fourth quadrant) , 
so that one cannot isolate the reversal there using the raw 
data of Figs |^a,d in the (RM, I) plane (Simard-Normandin & 
Kronberg 1979, 1980). Wavelet transformation provides ef- 
ficient scale filtering and thereby allows us to isolate weaker 
localized details seen against the background of an extended 
and stronger signal, and so the reversed field can be easily 
identified at both < I < 90° (the Sagittarius arm) and 
270° < I < 360° (the Carina arm). 

We note that the model with reversal presented in Fig. ^ 
reproduces well RM in the Cetus arc region. This is consis- 
tent with the fact that Loop II cannot be distinguished at 
low Galactic latitudes (Berkhuijsen 1971). 

Our fits to the data of Simard-Normandin et al. have 
a purely azimuthal magnetic field. We stress that this does 
not imply that the magnetic field in the Milky Way is purely 
azimuthal in reality, but merely indicates that its radial com- 
ponent cannot be estimated reliably from these data. There- 
fore, we prefer here a model with the minimum number of 
parameters; the data available do not warrant using more 
complicated (and realistic) models. For the same reason we 
do not employ a model of electron density including spiral 
arms (Taylor & Co rdes 1993). The data of Broten et al. 
(1988) used in Sect. 3.2.2 are not compatible with a purely 
azimuthal field, so we apply a more complicated model and 
estimate the pitch angle of magnetic field there. 



6.2.2 The catalogue of Broten et al. 

The catalogue of Broten et al. (1988) contains 69 low- 
latitude sources (21 more than the catalogue of Simard- 
Normandin et al.). There are 7 low-latitude sources in 
Simard-Normandin et al. which are absent in Broten et al. 
and which have also been included in the sample used in 
this section. A one-dimensional wavelet transform of these 
data (76 sources shown in Fig. 9a) is shown in Fig. 8a in the 
same format as in Fig. 0a. The large-scale structures due 
to the Orion arm, centred at I as 90° and 280°, are similar 
to those obtained from the data of Simard-Normandin et 
al. This is true also for the reversed magnetic field in the 
Sagittarius-Carina arm that produces structures at a scale 
35°-45°, slightly larger than from the high-latitude sources. 

A remarkable feature of Fig. 8a is that it has a new 
pair of structures centred at a scale less than 30°. This is a 
positive extremum at I « 125° seen against a negative back- 
ground produced by the local magnetic field in the Orion 
arm and a negative extremum at I ~ 190°. If these struc- 
tures can be attributed to the Perseus arm, magnetic field 
in that arm must be directed oppositely to the local field, 
in agreement with the earlier determination of Agafonov et 
al. (1988). The Perseus arm is seen at I = 145° -160° in syn- 
chrotron intensity and at I = 90°-140° in the optical range 
(Berkhuijsen 1971). The negative extremum attributed to 
the Perseus arm occurs close to the anticentre direction be- 
cause this relatively weak feature is best visible in the di- 
rection where the local contribution from the Orion arm is 
minimum. (Magnetic field reversal in the outer Galaxy was 
also advocated by Rand & Kulkarni 1989 and Clegg et al. 
1992.) 

The sample of 76 low-latitude sources includes some 
with very large values of |RM| (with the maximum of 
800radm~ 2 ). To assess the robustness of our results we 
removed sources with the largest |RM| reducing gradually 
the cut-off value. We show in Figs 9b and 8b the the orig- 
inal data and their wavelet transform for 68 sources with 
|RM| < 300radm -2 . The newly detected pair of structures 
becomes less pronounced, but remains well visible. 

The simplest model which is able to produce a similar 
pair of structures contains an additional reversal of magnetic 
field in the ou ter G alaxy. We first tried purely azimuthal 
field as in Sect. 3.2.1 but this failed to produce the required 



structures of the desired strength and location for any real- 
istic field strength. 

Therefore, we had to consider a more complicated model 
where magnetic lines are trailing logarithmic spirals. The 
simplest model of this kind includes three logarithmic spi- 
ral branches (for the Sagittarius-Carina, Orion and Perseus 
arms) with the pitch angle p = 15° and with magnetic field 
reversals between (or within) each arm. The value of the 
pitch angle adopted is close to 14° as inferred for the Galac- 
tic spiral arms from synchrotron intensity observations by 
Beuermann et al. (1985). With this model, we adjusted the 
strength of the magnetic field in each arm and the positions 
of the reversals. The best fit between the wavelet transforms 
for the original and model RM is obtained for the parame- 
ters given in Table 1 (for the Orion arm, B » 1.5 /iG was 
obtained) . The wavelet transform of the model RM is shown 
in Fig. 8c; we also show both the original and model RMs in 
Fig. 9 in the same format as in Fig. pi The structures due to 
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Figure 8. One-dimensional wavelet transforms for the combined 
set of 76 low-latitude sources, |b| < 10°, from Broten et al. (1988) 
and Simard-Normandin et al. (1981): (a) all the sources; (b) those 
with |RM| < 300°; and (c) model with two reversals and n c = ni 
as described in the text. Contour labels are in radm" 2 . 
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Figure 9. The RM data whose wavelet transforms are shown in 
the respective panels of Fig. 8. Curves shown in (c) have the same 
meaning as in Fig. pi 



the Perseus arm occur at 90° < I < 200°, where the sample 
of Broten et al. (1988) contains 13 additional sources. The 
errors in W and B given in Table 1 characterize the scatter 
of the values obtained from fitting to the data of Simard- 
Normandin et al. (Fig. 0a) and Broten et al. (Figs 8a,b). 

The histogram of ARM, the difference between the ob- 
served and model values of RM for the best fit model, is 
shown in Fig. [loj. The histogram is as close to a Gaussian as 
can be expected with 76 data points; the tail at positive val- 
ues is due to sources with |RM| > 300 radm -2 . The mean 
value of ARM in Fig. |l^ is 40 radm -2 for all 76 sources, 
but just 4radm -2 for sources with |RM| < 300radm -2 . 
The standard deviation of ARM is a ~ 100 radm -2 for 
the sources with |RM| < 300 radm -2 (and 180 radm -2 il 
the eight sources with |RM| > 300 radm -2 are included). 
A similar amount of scatter is implied by the observed 
rotation measure structure function, ranging from about 
40 radm -2 at a high latitude b = -20° (Minter & Spangler 
1996) to 200 radm -2 at low latitudes (0 < b < 4°) (Lazio, 
Spangler & Cordes 1990) at a scale 30°. A rough estimate 
of the scatter expected to arise within the Milky Way is 
a ~ Kn c B[{q 2 B +ql)dL] 1/2 ~ 40 rad m -2 (g| + g 2 ) 1/2 , where 
qb = 8B/B and q n = Sn e /n e are the relative fluctuations 




-300 -120 60 240 420 600 780 
ARM (rad m 2 ) 

Figure 10. The histogram of the difference between the observed 
RM of the low-latitude sources of Broten et al. (1988) and Simard- 
Normandin et al. (1981) and the best-field model RM with two 
magnetic field reversals; the bin width is 60 radm . 
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in magnetic field and electron density, d ~ 100 pc is the tur- 
bulent scale and L ~ 3 kpc is the pathlength (see Minter & 
Spangler 1996 for details). 

We note that the above fits including the Perseus arm 
place the reversal in the inner Galaxy at r = 6.5 kpc, i.e., 



somewhat farther from the Sun than the model of Sect. 3.2.1 



We consider the model of Sect. 5.2.1 to be more accurate in 
this respect and the corresponding estimate of the reversal 
radius is given in Table 1. 

The Orion arm can be different from the other arms 
being a short spur rather than a long arm; then its pitch 
angle can be significantly different from that of the Sagit- 
tarius and Perseus arms. If this is the case, the fit of this 
section will be distorted by our equal treatment of the three 
arms. Further refinement of the model will be useful when 
more extensive RM catalogues will become available. The 
presents fits look good enough without such refinements, so 
more complicated models cannot be justified. 

We do not see any signatures of the Scutum arm, an in- 
ner spiral arm next to the Sagittarius arm. Its longitudinal 
range is 24° I Z*> 43° (e.g., Beuermann et al. 1985, Vallee et 
al. 1988). Assuming the scale height of the regular magnetic 
field to be 1 kpc and the arm galactocentric distance of 5 kpc 
at I = 0, the distance to the tangential point is more than 
6.3 kpc; thus, the Scutum arm spans the latitude range of 
|6| ^ 9°. This window contains just two sources in the cata- 
logue of Broten et al. (1988), one with RM = 50radm -2 and 
the other with RM = — 840radm~ 2 , with the latter value of 
RM being too large to be attributed to the Galaxy. Vallee 
et al. (1988) claim that the Scutum arm produces a neg- 
ative contribution to RM using 9 sources 7 of which have 
negative RM. However, 6 of their sources with RM < 
have |6| > 10.4°. Moreover, two of them have RM = -699 
and — 840radm~ 2 ; such values of RM cannot be produced 
in the Milky Way. Vallee et al. justify their wider window 
in latitude, \b\ < 14°, by taking an overestimated magnetic 
field scale height of 1.25 kpc and underestimated distance to 
the tangent point (5 kpc). The scale height of the regular 
magnetic field is about lkpc in the disc near the Sun (e.g., 
Fletcher & Shukurov 2000) and even smaller in the inner 
Galaxy (Phillipps et al. 1981, Beuermann et al. 1985). 

6.3 The vertical symmetry of the large-scale 
magnetic field 

Although the structures in the RM distribution at the 
largest scale are slightly shifted to negative Galactic lati- 
tudes, they clearly indicate that the large-scale magnetic 
field in the Orion arm has the same direction above and be- 
low the Galactic midplane (see Sect. ^). This result agrees 
with most earlier studies of the Galactic magnetic field (with 
a few exceptions - see Sect. ^). However, earlier results were 
based either on a naked-eye inspection of noisy RM maps 
or on fitting magnetic field models having presupposed even 
symmetry. Our results provide a rare objective confirmation 
of the anticipated even symmetry of the horizontal magnetic 
field in the Orion arm of the Milky Way. 

The fact that the field direction is the same above and 
below the Galactic equator is compatible with the overall 
quadrupole symmetry of magnetic fields generated by disc 
dynamos. However, this cannot be considered as a com- 
pletely convincing proof of the dynamo origin of the Galac- 



tic magnetic field because the dynamo-generated quadrupole 
magnetic field must also have a vertical magnetic field an- 
tisymmetric with respect to the midplane. Hence, strictly 
speaking, one has to measure B z before it can be con- 
cluded definitely that the observed symmetry directly con- 
firms galactic dynamo models. Still, the even symmetry of 
the azimuthal field proved here significantly restricts theo- 
ries of magnetic field origin (see Beck et al. 1996). 

Even when the contribution of Loop I has been sub- 
tracted, the centres of the dominant |RM| maxima for high- 
latitude sources remain shifted from the midplane by about 
15° to the south. (We have verified that the overall sym- 
metry remains even when we discard also the data in a re- 
gion in the Southern hemisphere located symmetrically to 
Loop I.) This residual asymmetry can hardly be explained 
by the shift of the Sun's position away from the Galactic 
midplane. A displacement of the Sun by about 400 pc to the 
north would be required to reproduce the observed residual 
asymmetry, whereas the actual displacement is as small as 
50 pc. 

Another cause of the asymmetry could be a local con- 
tamination from Loop II, centred at (l,b) = (110°,— 32°), 
whose radius is 50° (Berkhuijsen 1971). However, there are 
no local magneto-ionic structures of a similar size located 
symmetrically to Loop II in the fourth quadrant (the Gum 
nebula has a small angular radius of 18° and occurs at b = 
- e.g., Vallee 1993). Therefore, the symmetry of the residual 
shift in the large-scale wavelet extrema between the first and 
fourth Galactic quadrants does not support this idea. 

The residual shift of the centres of the dominant RM 
structures to negative Galactic latitudes can be a result of 
a genuine global vertical asymmetry of the magnetic field 
with a stronger large-scale field in the southern hemisphere. 
The possibility of an asymmetry of this type was suggested 
by Sokoloff & Shukurov (1990) as a consequence of the odd 
symmetry of magnetic field in the Galactic halo (see also 
Brandenburg et al. 1992; Brandenburg, Moss & Shukurov 
1995). The odd symmetry of the halo magnetic field has been 
suggested by the spherical harmonic analysis of extragalactic 
RM by Seymour (1984). The shift to b ~ -15° then implies 
that the vertical distribution of the field has a maximum 
at z = — Lstanbo ~ 400 pc. We stress that this estimate 
is only illustrative and has an uncertainty of at least 50 
percent. 

This interpretation is not unique. A similar asymmetry 
would arise due to a local perturbation in magnetic field 
and/or electron density. We can estimate the implied prop- 
erties of such a perturbation assuming, for example, that the 
Sun is located near the top of a magnetic loop extending ver- 
tically above the midplane. Such a loop can result from the 
Parker instability. If the base of the loop is at the midplane, 
2 = 0, and the Sun is at the height zq — 50 pc, then the 
horizontal size of the loop is L ~ 2z0/|tan&o| — 400 pc. 
The difference between the values of w at b = —15° and 
b = +15° for I = 98° in Fig. |(i is Aw ~ 4radm -2 . Then 
the field enhancement in the loop (relative to the average 
field) is estimated as AB ~ Aw/(Kn c L) ~ 0.5 fj,G. 
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7 CONCLUSIONS 

Using the technique of wavelet analysis adapted to nonuni- 
form data grids, we have studied the large-scale Galactic 
magnetic field using Faraday rotation measures of extra- 
galactic radio sources. Our main results are based on the 
catalogues of Simard-Normandin et al. (1981) and Broten 
et al. (1988). We have also incorporated some measurements 
from other catalogues into a joint catalogue containing 841 
source. 

We have tested the stability of the results with respect 
to modifications of the data sample. In particular, results 
obtained from the joint catalogue are rather unstable con- 
firming that the uniformity of the data sample is important. 
We use other catalogues to select unreliable measurements in 
the catalogue of Simard-Normandin et al. (1981) and discard 
only four sources (see Sect. |^) ; this agrees well with the value 
expected from statistical arguments (P.P. Kronberg, private 
communication). Therefore we used only 551 sources from 
the catalogue of Simard-Normandin et al. (1981). 

The large-scale magnetic field in the local spiral arm is 
dominated by the azimuthal component. We confirm that 
there is a reversal in the magnetic field in the inner Galaxy 
at a distance of 0.6-1 kpc from the Sun. Using an axisym- 
metric model of the Galactic distribution of thermal elec- 
trons (Cordes et al. 1991), we estimate the strength of the 
magnetic field in the Orion arm as 1.4 ± 0.3 /iG. 

The extension of the magnetic fields reversal between 
the Orion and Sagittarius arms to the Carina arm region 
at negative Galactic latitudes has been detected. Faraday 
rotation measures of low-latitude sources from the catalogue 
of Broten et al. (1988) are consistent with a field reversal 
between the Orion and Perseus arms. The average local pitch 
angle of the magnetic field is about 15°. 

The extent, in the azimuthal direction, of the region 
with reversed large-scale magnetic field in the Solar vicin- 
ity cannot be determined from the available data. We can 
only set its lower limit at 3 kpc. Therefore, the region can 
be relatively small not only in radius (2-3 kpc) but also in 
azimuth, rather than extended over the whole circle around 
the Galactic centre (if the field is axisymmetric) or along a 
spiral line (for a bisymmetric field) (Shukurov 2000). Mor- 
phologically, this would be in accord with the Orion arm be- 
ing a short spur between the Sagittarius and Perseus arms - 
the field would be reversed in this spur. A similar magnetic 
structure has been found in the galaxy M51 by Berkhuijsen 
et al. (1997), where the regular magnetic fields in the disk 
is reversed in a region of a size about 8 kpc in azimuth and 
3 kpc in radius. 

The dominant even symmetry of the large-scale local 
magnetic field with respect to the Galactic equator has been 
demonstrated convincingly using wavelets. This symmetry 
becomes obvious after filtering out RM structures at angular 
scales smaller than about 70° (these smaller-scale structures 
include Loop I). 

Our results are consistent with a moderate asymmetry 
of the local large-scale Galactic magnetic field with respect 
to the equator (with the |RM| maxima occurring at about 
b — —15°) that can be explained by stronger field in the 
southern Galactic hemisphere or by the Sun being located 
close to the top of a magnetic loop such as that arising from 



the Parker instability. The horizontal size of the asymmetric 
region is at least 3 kpc if it is global and 400 pc if it is local. 

Loop I affects the values of RM within a region which 
extends to I ~ 50°, that is by 10°-20° farther than the edge 
of Loop I as seen in synchrotron emission. 



APPENDIX A: WAVELET ANALYSIS FOR 
DISCRETE DATA IRREGULARLY 
DISTRIBUTED ON A SPHERE 

A basic property of wavelets is that their average value van- 
ishes, see Eq. (fa). In an ideal situation and for a flat space 
this requirement is trivially ensured by a proper choice of 
the functional form of the analyzing wavelet il>(x). However, 
a curved space involves an intrinsic parameter having the di- 
mension of length, namely the local curvature radius. There- 
fore, it is impossible to introduce an exactly self-similar fam- 
ily of wavelets in a curved space, for example on a sphere. 
The intrinsic deviations from self-similarity are significant 
when the wavelet scale a is comparable to or larger than 
the curvature radius. Here we suggest a modification of the 
wavelet transformation technique to ensure an approximate 
self-similarity which becomes almost exact when the scale 
of the wavelet a is much smaller than the curvature radius. 

Another problem is that one should take special care 
about condition (|^) in numerical realizations of the wavelet 
transformation, especially when the data points are irregu- 
larly distributed in space. It is nontrivial to ensure that the 
average value of the wavelet of a scale a is zero when a is 
comparable to or smaller than the largest separation of the 
data points. 

The two difficulties mentioned above can be resolved 
simultaneously using the gapped (or adaptive) wavelet tech- 
nique suggested by Frick et al. (1997a) for a one-dimensional 
case (see also Frick, Grossmann & Tchamichian 1998). Here 
we generalize this technique to wavelets defined on a spher- 
ical surface. 

Let us represent the analyzing wavelet in the form 

ip(x) = h{x)<f>{x) , (Al) 

where &(x) is a positive-definite envelope and h(x) repre- 
sents the filling of the wavelet. In the case of the 'Mexi- 
can hat' (^) in a plane we have $(a;) = exp (— ^a; 2 ) and 
h(x) = 2-x 2 . 

Following Frick et al. (1997a), we use the wavelet of the 
form 

ifj a , r (r) = a" 1 [h(s/a) - K{a, r)] $(s/a) (A2) 

in the numerical realization of the wavelet transformation. 
Here s, the angular separation of r and r', is a function of 
r = (I, b), the position of the centre of the wavelet and r' = 
(l',b'), a current position on the unit sphere (r and r' are 
defined on the unit sphere and should not be confused with 
radius vectors). The quantity )C(a, r) is determined from the 
condition }] . ipi = 0, where ipi are the values of the wavelet 
at the data grid points (li,bi). This yields 

K(a,l,b) = \^2$(si/a)\ ^2h(si/a)$(si/a) , (A3) 

i i 

where Si are the distances of the data points to the centre of 
the wavelet (/, b). Thus, we have c = 2 — PC(a, r) in Eq. (^). 
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Figure Al. R as a function of a for a sphere (dashed line) and 
for a plane (solid). Error bars indicate discretization errors for 
the catalogue of Simard-Normandin et al. (1981). 



All the integrals (g) were calculated on the irregular 
data grid determined by the particular catalogue using the 
Monte-Carlo technique as 



w(a, I, b) = AS V'o.r 0;) RM 



(A4) 



where summation is carried out over the data points. Here 
-l 



AS 



is the area element. 



A significant dependence of IC(a, r) on position for small 
a affects the scale resolution of the wavelet technique. In 
terms of the radius R at which a wavelet has the first zero, 
R = o[2 — K,(a, r)] 1 / 2 , the resulting wavelet transform is un- 
reliable at those scales where the scatter in R arising from 
the position variations of K(a, r) is comparable with a. The 
specific form of )C(a, r) depends on the distribution of the 
data poi nts on the sphere, i.e. on the catalogue chosen. In 
Fig. 



Al 



we show R as a function of a for the catalogue 
of Simard-Normandin et al. (1981) and indicate the scat- 
ter arising from the inhomogeneous distribution of the data 
points. For R 20°, the relative error exceeds 25%; we 
therefore consider the results obtained from this catalogue 
to be reliable at scales R > 20°. The catalogue of Broten 
et al. (1988) contains more sources than that of Simard- 
Normandin et al. (1981). However, the sources in the former 
catalogue are distributed less uniformly in the sky, so the re- 
sulting discretization error is even larger than for the latter 
catalogue. Therefore, the resolution limit of both catalogues 
is adopted as R = 20°. Note that this restriction does not 



apply in Sect. 6.2 where we consider low-latitude sources 



and the separation of data points is reduced by projection 
on the equator. 

We have tested our implementation of the wavelet trans- 
formation using artificial data sets given on both a regular 
grid in (l,b) (with a 81 x 41 mesh) and the grids of the 
catalogues considered in this paper. We specified a test dis- 
tribution of RM in the sky consisting of four structures at 
different scales: two spots of different signs having a scale 
R = 70° and located at (I, b) = (±90°, 0), one positive spot 
of a scale R = 40° at (l,b) = (45°, -45°), and another 
negative spot of R — 25° located at the Galactic North- 
ern pole. The amplitudes of all the structures were equal to 



each other. In addition, we superimposed a white noise of a 
relative amplitude 15%. 

We calculated the values of RM of the test distribution 
at those positions in the sky where Simard-Normandin et al. 
(1981) give their values of RM and also on the grid of the 
joint catalogue. We then performed wavelet transformation 
of the resulting discrete data sets specified on the irregular 
grids and checked that the wavelet transform successfully 
isolates the structures in the input signal, and that these 
structures, if well separated in scale, result in local maxima 
of the integral energy spectrum M(a) as given by Eq. (m) 
(calculated for the same discrete data set). 

This test confirms that our algorithm correctly recovers 
structures on real data grids: all the four structures can be 
clearly seen in the wavelet transform maps. Traces of the 
structures at larger scales can also be seen at smaller scales 
because our wavelet provides a good spatial resolution (i.e., 
it is sensitive to the position of a structure) and only a mod- 
erate scale resolution. The results remain stable when white 
noise has been added to the data. 
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